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Abstract 

The transverse target spin azimuthal asymmetry j4^"y ™in hard exclusive production of p° me- 
sons was measured at COMPASS by scattering 160GeV/c muons off transversely polarised protons 
and deuterons. The measured asymmetry is sensitive to the nucleon helicity-flip generalised parton 
distributions E q , which are related to the orbital angular momentum of quarks in the nucleon. The 
Q 2 , xgj and pj, dependence of A s ^f ^ s hs presented in a wide kinematic range: 1 (GeV/c) 2 < 
Q 2 < 10(GeV/c) 2 , 0.003 < x Bj < 0.3 and 0.05(GeV/c) 2 < p% < 0.5(GeV/c) 2 for protons or 
0.10(GeV/c) 2 < p^r < 0.5(GeV/c) 2 for deuterons. Results for deuterons are obtained for the 
first time. The measured asymmetry is small in the whole kinematic range for both protons and 
deuterons, which is consistent with the theoretical interpretation that contributions from GPDs E" 
and E d approximately cancel. 
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1 INTRODUCTION 



1 Introduction 

Hard exclusive electro- and muoproduction of mesons on nucleons has played an important role in studies 
of strong interactions and recently gained renewed interest as it allows access to generalised parton 
distributions (GPDs) 0343]]. The GPDs provide a novel and comprehensive description of the partonic 
structure of the nucleon and contain a wealth of new information. In particular, they embody both 
nucleon electromagnetic form factors and parton distribution functions. Furthermore, GPDs provide a 
description of the nucleon as an extended object, referred to as 3-dimensional nucleon tomography [6]- 
8], which correlates longitudinal momenta and transverse spatial degrees of freedom of partons. The 
evaluation of GPDs may for the first time provide insight into angular momenta of quarks, another 
fundamental property of the nucleon ElO. The mapping of nucleon GPDs, which very recently became 
one of the key objectives of hadron physics, requires a comprehensive program of measuring various hard 
exclusive processes in a broad kinematic range, in particular deeply virtual Compton scattering (DVCS). 
Hard exclusive meson production provides independent and complementary information. 

In perturbative QCD (pQCD), there exists a general proof of factorisation [ 9 ] for exclusive meson produc- 
tion by longitudinal virtual photons. In this case the amplitude for hard exclusive meson leptoproduction 
can be factorised into a hard-scattering part and soft parts, the latter depending on the structure of the 
nucleon described by GPDs and on the structure of the meson described by its distribution amplitude 
(DA). No similar proof of factorisation exists for transverse virtual photons. However, pQCD-inspired 
models taking into account parton transverse momenta have been proposed lflQl4T2l . which describe rea- 
sonably well the behaviour of the cross sections for both longitudinal and transverse photons, <tl and cjt, 
respectively. 

At leading twist, meson production is described by four types of GPDs: H* , E*, H*, E' , where / de- 
notes a quark of a given flavour or a gluon. The GPDs are functions of t, x and £, where t is the squared 
four-momentum transfer to the nucleon, x the average and £ half the difference of the longitudinal mo- 
menta carried by the struck parton in the initial and final states. In addition, there is a scale dependence 
of GPDs which is not explicitly shown here. Depending on the quark content and the quantum numbers 
of the meson, there exists sensitivity to various types of GPDs and different quark flavours. In particular, 
production of vector mesons is sensitive only to GPDs H* and E* . 

The GPDs attracted much attention after it was shown that the total angular momentum of a given parton 
species / is related to the second moment of the sum of GPDs H* and E$ via the Ji relation O : 



which holds for any value of £. The spin-independent cross sections for DVCS and for vector meson pro- 
duction by longitudinal photons on a proton target are mostly sensitive to the nucleon-helicity-conserving 
GPDs Hf, with GPDs being suppressed in the COMPASS kinematic domain. However, the GPDs 
Ef are of special interest, as they are related to the orbital angular momentum of quarks. They describe 
transitions with nucleon helicity flip, in which orbital angular momentum must be involved due to total 
angular momentum conservation. It was pointed out that the spin dependent cross sections for DVCS on 
transversely polarised protons lfT3l [141 and for exclusive vector meson production by longitudinal pho- 
tons on transversely polarised nucleons lfl"5ll are sensitive to the 'elusive' nucleon helicity-flip GPDs E* . 
Access to GPDs E? is also possible by measurements of the cross section for DVCS on an unpolarised 
neutron target lfl4l . 

Measurements of the lepton helicity dependent DVCS cross section on neutrons were performed by the 
JLAB Hall A collaboration [16] and the transverse target spin asymmetries for DVCS from transversely 
polarised protons were measured by the HERMES experiment [17]. Model-dependent estimates of the 
total angular momenta of quarks, J u and J d , derived from the results of these measurements, indicate a 
large value for J u and a value close to zero for J d , in agreement with results from lattice QCD [181- 




(1) 
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In exclusive vector meson production on transversely polarised targets the observable sensitive to the 
GPDs E$ is the azimuthal asymmetry A^j? ^ s \see Sec. [2] for the definition). Here, the indices U and 
T refer to the beam spin independent and transverse target spin dependent cross section, and sin(^> — <pg) 
indicates the type of azimuthal modulation of the cross section. The GPDs E$ appear at leading twist 
only in this azimuthal asymmetry for vector meson production by longitudinal photons. 

The only previous measurement of A^^^for exclusive p° electroproduction on transversely polarised 
protons was performed by the HERMES experiment ffT9ll20l . Its separate extraction for longitudinally 
and transversely polarised p° mesons gave values consistent with zero. A model-dependent attempt was 
made [19] to extract the value of the total angular momentum J u of u quarks in the proton, although 
limited by large experimental uncertainties. 

In this paper, we present results on A s ^f ^for exclusive p° meson muoproduction on transversely po- 
larised protons and deuterons. The experiment was carried out at CERN by the COMPASS collaboration 
using the 160GeV/c muon beam and a polarised target filled either with lithium deuteride ( 6 LiD) or 
ammonia (NH3) to provide polarised deuterons or protons, respectively. 



2 Theoretical framework 

The cross section of hard exclusive p° leptoproduction, p N — > p p° N, on a transversely polarised nu- 
cleon depends on the photon virtuality Q 2 , the Bjorken variable XBj, t, (ft and 4>s ED - Here <p is the 
azimuthal angle between the lepton scattering plane and the plane containing the virtual photon and the 
produced meson (hadron plane), while <ps is the azimuthal angle of the target spin vector around the 
virtual photon direction relative to the lepton scattering plane (see Fig. [I]). A summary of the kinematic 




Fig. 1: Kinematics of exclusive meson production in the target rest frame. Here k, k', q and v represent 
three-momentum vectors of the incident and the scattered muons, the virtual photon and the meson. St 
is the component of the target spin vector S (not shown) perpendicular to the virtual photon direction. 



variables used in this paper is given in Table [T] 

In the COMPASS kinematic region the cross section can be expressed as: 

? — 1 
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where only terms relevant for the present analysis are shown explicitly. The general formula for the cross 
section for meson leptoproduction, which contains the dependence on the projectile and target spins and 
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2 THEORETICAL FRAMEWORK 



Table 1: Kinematic variables. 



k four-momentum of incident muon 

k' four-momentum of scattered muon 

p four-momentum of target nucleon 

v four-momentum of p° meson 

q — k — k' four-momentum of virtual photon 

Q 2 = — q 2 invariant negative mass squared of virtual photon 

W = \J (p + q) 2 invariant mass of the 7* — N system 

M p proton mass 

v = (p-q) /M p energy of virtual photon in the laboratory system 

xjjj = Q 2 J (2M p v) Bjorken scaling variable 

y = (p-q)/(p- k) fraction of lepton energy lost in the laboratory system 

M n „ invariant mass of tt + tt~ system 

t=(q — v) 2 square of the four-momentum transfer to the target nucleon 

pir transverse momentum squared of vector meson with 

respect to the virtual photon direction 

E p o energy of p° in the laboratory system 

M x = (p + q — v) 2 missing mass squared of the undetected system 



-Emiss = {{p + q — v ) 2 —p 2 )/(2Mp) missing energy of the undetected system 

= is-E p o+t/(2M p ) 



the complete azimuthal dependence, can be found in Ref. [21 ]. The component of the transverse target 
spin perpendicular to the virtual photon direction, St, is in the COMPASS kinematic region very well 
approximated by the corresponding component perpendicular to the direction of the incoming muon. 
The virtual photon polarisation parameter e is given by 



i-y-?y 2 7 2 



where y is the virtual-photon fractional energy, 7 = 2xBjM p /Q and M p the proton mass. 



(3) 



The symbols o % mn represent spin-dependent photoabsorption cross sections or interference terms, which 
are proportional to bilinear combinations of helicity amplitudes A l m for the subprocess 7* iV — > p°N, 



oc 



E 

spins 



(•An) -^ni (4) 



where the dependence on kinematic variables is omitted for brevity. The virtual-photon helicity is de- 
noted by m(n) = 0, ±1, the target nucleon helicity by i(j) = ±5, and the notation is restricted to 0, +, — 
for legibility. 

The transverse target spin dependent cross section is accessed experimentally by measuring the azimuthal 
asymmetry A S ™ T ^ which is proportional to the sin(^> — ^s) moment of the cross section for a trans- 
versely polarised target. It can be expressed as 

A sm(<p-^ s ) = Im(g++ + e Pro" ) = Im K; + e gaT ) (5) 

The denominator contains the spin-averaged cross section with contributions from both transverse and 
longitudinal virtual photons. The leading twist interference term corresponds to p° production by 
longitudinal photons, while the higher twist term <r++ corresponds to production by transverse photons. 
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The former term is proportional to a weighted sum of convolutions of the GPDs E$ with the DA of the 
produced meson and a hard scattering kernel [ 15"]. The weights depend on the contributions of quarks of 
various flavours and gluons to the production of a given vector meson. 

The direct method to separate terms arising from production by longitudinal and transverse photons 
is the Rosenbluth separation. However, it is not feasible with the present data as only measurements 
at one beam energy are available. In Ref. ll22l another method was proposed which can be used for 
vector meson production in the approximation of s-channel helicity conservation and exploiting the 
decay angular distributions of the meson. It was applied by HERMES for the analysis of p° production 
on transversely polarised protons. In our analysis we do not attempt such a separation. Nevertheless, 
the present 'unseparated' results can be compared to the predictions of existing models which take into 
account also higher twist effects. 



3 Experimental set-up 

The experiment E3ll was performed using the high intensity positive muon beam from the CERN SPS 
M2 beam line. The instantaneous /i + beam intensity during extraction is about 4 • 10 7 /s. The average 
beam momentum is 160GeV/c with a spread of 5GeV/c. The momentum of each incident muon is 
measured upstream of the experimental area with a relative precision better than 1 %. The p + beam is 
longitudinally polarised by the weak decays of the parent hadrons. Note that the beam polarisation does 
not affect the measurement of A s ™j? ^ • 

The beam traverses a polarised solid-state target which contains 120 cm total length of polarisable ma- 
terial, which is either NH3 for polarised protons or 6 LiD for polarised deuterons. Both protons and 
deuterons can be polarised either longitudinally or transversely with respect to the beam direction. A 
mixture of liquid 3 He and 4 He, used to refrigerate the target, and a small admixture of other nuclei are 
also present in the target. It consists of either three (NH3) or two ( 6 LiD) separate cells with polarisable 
material, placed one after another along the beam. The spin directions in neighbouring cells are op- 
posite. Both target configurations allow for a simultaneous measurement of azimuthal asymmetries for 
the two target spin directions to compensate flux-dependent systematic uncertainties. In order to reduce 
systematic effects of the acceptance, the spin directions are reversed periodically about every week with 
polarisation measurements before and after reversal. The three-cell configuration results in a more bal- 
anced acceptance for cells with opposite polarisation, which leads to a further reduction of systematic 
effects. The achieved polarisation, Pt, is about 0.8 for protons (NH3) and 0.5 for deuterons ( 6 LiD) with 
relative uncertainties of 3 % and 5 %, respectively. 

The fraction of polarisable material in the target weighted by the corresponding cross sections is quanti- 
fied by the dilution factor, /, that depends on the considered reaction. It is calculated using the measured 
material composition of the target and the nuclear dependence of the cross section for the studied reac- 
tion. For incoherent exclusive p° production the dilution factor is typically 0.25 for the NH3 target and 
0.45 for the 6 LiD target. See Sec.|6]for more details. 

The target is housed in a large superconducting solenoid providing a field of 2.5 T along the beam direc- 
tion. From 2002 to 2004 the angular acceptance was ±70 mrad at the upstream edge of the target. From 
2006 onwards an upgraded target magnet with a new large-aperture solenoid was used. It provides an 
angular acceptance of ±180 mrad for the upstream target edge resulting in an increased hadron accep- 
tance. The transverse holding field of up to 0.5 T, provided by a dipole coil, is used for adiabatic spin 
rotation and for measurements with transverse target polarisation. 

The COMPASS spectrometer is designed to reconstruct scattered muons and produced hadrons in wide 
momentum and angular ranges. It consists of two stages, each equipped with a dipole magnet, to measure 
tracks with large and small momenta, respectively. In the high-flux region, in or close to the beam, 
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4 EVENT SAMPLE 



tracking is provided by stations of scintillating fibres, silicon detectors, micromesh gaseous chambers 
and gas electron multiplier chambers. Large-angle tracking devices are multiwire proportional chambers, 
drift chambers and straw detectors. Muons are identified in large-area mini drift tubes and drift tubes 
placed downstream of hadron absorbers. Each stage of the spectrometer contains an electromagnetic and 
a hadron calorimeter. The identification of charged particles is possible with a RICH detector, although 
in this analysis we have not utilised the information from the RICH. 

The data recording system is activated by several triggers. For inclusive triggers, the scattered muon is 
identified by a coincidence of signals from trigger hodoscopes. Semi-inclusive triggers select events with 
a scattered muon and an energy deposit in a hadron calorimeter exceeding a given threshold. Moreover, a 
pure calorimeter trigger with a high energy threshold was implemented to extend the acceptance towards 
high Q 2 and xsj- In order to suppress triggers due to halo muons, veto counters upstream of the target 
are used. The COMPASS trigger system covers a wide range of Q 2 , from quasi-real photoproduction to 
deep inelastic interactions. 

4 Event sample 

The results presented in this paper are based on the data taken with the transversely polarised 6 LiD target 
in 2003-2004 and with the transversely polarised NH3 target in 2007 and 2010. The phase space of 
the incoming beam is equalised for all target cells using appropriate cuts on position and angle of beam 
tracks. An event to be accepted for further analysis is required to have an incident muon track, a scattered 
muon track and exactly two additional tracks of oppositely charged hadrons, all associated to a vertex 
in the polarised target material. Figure [2] shows the distribution of the reconstructed vertex position zy 
along the beam axis. In this figure as well as in Figs. [3] to [5] the distributions are obtained applying all 
cuts except those corresponding to the displayed variable. 




Fig. 2: Distributions of the z-coordinate of the primary vertex zy for NH3 (left) and 6 LiD data (right). 
The vertical lines indicate the applied zy cuts. 



In order to obtain a data sample in the deep inelastic scattering region, the following kinematic cuts are 
applied: 1 (GeV/c) 2 < Q 2 < 10(GeV/c) 2 , where the upper limit is chosen to remove the region of Q 2 
where the fraction of non-exclusive background is large; 0.1 < y < 0.9, in order to remove events with 
large radiative corrections (large y) or poorly reconstructed kinematics (low y). The second cut removes 
also events from the region of hadron resonances at small values of W. A small residual number of such 
events is removed by requiring W to be larger than 5 GeV/c 2 . 
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As RICH information is not used in this analysis, the charged pion mass hypothesis is assigned to each 
hadron track. In order to select events of incoherent exclusive p° production the following additional 
cuts are applied, which will be justified below: 

0.5GeV/c 2 < M n7T < l.lGeV/c 2 , (6) 

- 2.5 GeV < E miss < 2.5 GeV, (7) 

E p0 > 15 GeV, (8) 

and 

0.1 (GeV/c) 2 <Pt< 0.5 (GeV/c) 2 for 6 LiD (9) 

or 

0.05(GeV/c) 2 <p^<0.5(GeV/c) 2 forNH 3 . (10) 

These cuts allow us to minimise the effects of various types of backgrounds such as: (i) semi-inclusive 
deep-inelastic (SIDIS) production of a p° meson or 7r + 7r _ pair, (ii) p° production with diffractive disso- 
ciation of the target nucleon, (hi) exclusive non-resonant tt + it~ pair production, (iv) coherent exclusive 
p° (or non-resonant 7r + 7r~ pair) production on a target nucleus. 

Figure [3] shows the distributions of M n7T for the NH3 and 6 LiD targets. A clear peak of the p° resonance 
is visible on top of a background arising from (i) and (hi). The selection on M n7T (Eq. ([6])) is optimised to 
minimise the effect of exclusive non-resonant tt + it~ pair production (hi) that will be discussed in Sec. [5] 




Fig. 3: Distributions of M vir for the NH3 (black, left scale) and 6 LiD (red, right scale) data. Vertical lines 
indicate the applied cuts. In order to exclude events originating from production of <fi mesons decaying 
into two charged kaons, the cut Mkk > 1.04GeV/c 2 is applied, where Mkk is the invariant mass of 
the two hadron system calculated assuming that both hadrons are kaons. 



Because slow particles are not detected, exclusive events are selected by the cut on missing energy given 
by Eq. ([7]). The selected range is referred to as 'signal region' in the following. In the -E m i ss distributions 
presented in Fig. [4] the peak at EWss ~ is the signal of exclusive p° production. The width (rms) of 
the peak, a rj 1.25 GeV, is due to spectrometer resolution which motivates the cut on E^ss (Eq. (|7])). 
Non-exclusive events (i) and (ii), where in addition to the recoil nucleon other undetected hadrons are 
produced, appear at EWss above about zero. However, due to the finite resolution they cannot be resolved 
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4 EVENT SAMPLE 
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Fig. 4: Distributions of -E m i ss for the NH3 (black, left scale) and 6 LiD (red, right scale) data. Vertical 
lines indicate the applied cuts to select the signal. 
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Fig. 5: Distributions of p\ for the NH3 (black, left scale) and 6 LiD (red, right scale) data. Vertical lines 
indicate the applied cuts. 



from the exclusive peak. The observed difference of -E m i ss distribution shapes between the two samples at 
large E'miss is due to the increase of the angular acceptance of the COMPASS setup mentioned in Sec. [3] 

The p\ distributions are shown in Fig. [5] We choose to use p\ rather than t or t' = \t\ — to, where to 
is the minimal kinematically allowed \t\, because in the COMPASS kinematic region px is determined 
with better precision by a factor of 2 to 5. In addition, the t' distribution is distorted because to, which 
depends on W, Q 2 , and M\, is poorly determined for non-exclusive background events [24]. The 
shown prr distributions indicate at small p\ values contributions from coherent p° production on target 
nuclei. Coherent events are suppressed by applying the lower cuts given by Eqs ( 9|l0 l. A study of p\ 
distributions shows that in addition to exclusive coherent and incoherent p° production a third component, 
originating from non-exclusive background (i), is also present and its contribution increases with p^, thus 
requiring also an upper cut. Therefore, in order to select the sample of events from incoherent exclusive 
p° production, the afore-mentioned p\, cuts were applied, which are indicated by vertical lines in Fig. [5] 
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After all selections the final samples for incoherent exclusive p production consist of about 797000 for 
the NH3 target and 97000 events for 6 LiD target. The mean values of the kinematic variables Q 2 , xbj, 
y, W and p\ are given in Table [2J 

Table 2: Mean values of the kinematic variables for proton and deuteron data. 





(Q 2 ) (GeV/c) 2 




(y) 


(W) (GeV/c 2 ) 


(p 2 T ) (GeV/c 2 ) 


proton data 


2.2 


0.039 


0.24 


8.1 


0.18 


deuteron data 


2.0 


0.032 


0.27 


8.6 


0.23 



5 Background to exclusive p° production 

The most important background contributions introduced in the previous section are discussed here in 
more detail. 

(i) The SIDIS contribution constitutes the largest source of background to the exclusive sample. It is 
estimated using Monte Carlo (MC) samples generated by the LEPTO generator with the COMPASS 
tuning of JETSET parameters l25l . The detector response is simulated using the description of either 
the NH3 or the 6 LiD set-up for transverse target polarisation. Simulated data are subject to the same 
selection criteria as real data and analysed in the same bins of kinematic variables. Comparing the E m i ss 
distributions of real and simulated data reveals insufficient agreement for E m iss > 7 GeV, where only 
SIDIS background contributes. The situation is improved considerably by weighting the h + h~ MC data 
in every E m iss bin i, by the ratio of numbers of like-sign events from real and MC data: 

W i N h+h+ + N h-h-> (Uj 

independent of other kinematic variables. As like-sign data contain only background, the weighting can 
be applied over the full Emiss range. This weighting procedure relies on the assumption that weights 
obtained from like-sign data are applicable to unlike-sign data. This assumption is supported by the 
observation that u>' lke ~ u>™ llke = N^£ al /N^ c holds at large E m i ss , despite of different shapes of 
Emiss distributions. The shape of the resulting weighted E m j ss distribution for unlike-sign MC data is 
parameterised for each individual target cell in every bin of Q 2 , %sj or p\. As the acceptance does not 
show a 4> — 4>s dependence, the MC data is not binned in this variable. 

For the determination of the asymmetries as described in Sec.|6] the NH3 and 6 LiD real data are binned 
in the same way in Q 2 , xsj or p\ per target cell, and also in (/) — (/) s and according to the target spin 
orientation. Here, a binning in 4> — 4>s is preferred over a simultaneous binning in 4> and <pg due to 
lack of statistics. In every such bin, the Emiss distribution is fitted using a Gaussian for the signal of 
exclusive events in conjunction with the above explained fixed shape for the SIDIS background, with 
free normalisation. The data are corrected for this background on a bin-by-bin basis. As an example, 
Fig. [6] illustrates the two-component fit for the Q 2 bin with the largest background contribution. From 
these fits, the fractions / s id; s of SIDIS events in the signal region (see Eq. ([7])) are found to vary between 
0.05 and 0.4, depending on kinematics and target set-up. On average, / s idi s equals 0.18 for deuteron and 
0.22 for proton data. 

(ii) Diffractive dissociation of the target nucleon into several particles is another type of background. 
Without a recoil detector, such events cannot be resolved from exclusive events by requirements on 
missing energy unless the mass Mx of the recoiling system exceeds 2.4GeV/c 2 . The high-Mx tail of 
diffractive dissociation events can be seen in Fig. [6] next to the exclusive peak as a small enhancement 
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5 BACKGROUND TO EXCLUSIVE p° PRODUCTION 




Fig. 6: The E miss distributions in the range 2.4(GeV/c) 2 < Q 2 < 10.0(GeV/c) 2 , together with signal 
plus background fits (solid curves) for the NH3 (left) and 6 LiD (right) samples. The dotted and dashed 
curves represent the signal and background contributions, respectively. 



over the SIDIS background. Diffractive dissociation background is examined using a MC event gen- 
erator called HEPGEN [26]. This generator is dedicated to studies of hard exclusive single photon or 
meson production processes in the COMPASS kinematic domain. In addition, it allows also to generate 
single photon or meson production accompanied by diffractive dissociation of the nucleon. Using both 
exclusive and nucleon-dissociative p° events generated by HEPGEN, which are reconstructed and se- 
lected as the real data, the contribution from low-mass diffractive dissociation of the nucleon is found to 
be rj 14% of the incoherent exclusive p° signal. No attempt is made to remove this type of background, 
motivated by HERA results on p° production where for unpolarised protons the angular distributions 
of proton-dissociative events were found consistent with those of exclusive events [ 27] - l29ll . Using the 
i^miss shape from HEPGEN, three-component fits to the experimental -Eniiss distributions show negligible 
impact of nucleon-dissociative events onto the determination of / s idi s - 

(iii) The background contribution due to exclusive non-resonant tt + 7t~ production, in particular the im- 
pact of its interference with resonant p° — > 7r + ir~ production, is studied in bins of Q 2 , XBj or p\, while 
suppressing the SIDIS background by a restrictive cut on missing energy, — 2.5 GeV < < OGeV. 
As an example, the M n+7T - distribution is shown in Fig. [7] for a selected Q 2 range, l(GeV/c) 2 < 
Q 2 < 1.2(GeV/c) 2 . The modification of the p° relativistic Breit-Wigner shape in the presence of non- 
resonant 7r + 7r~ events, which is observed for the distribution, is taken into account by applying either 
the Soding ll30ll or the Ross-Stodolsky PTll approach. The result of the fit to the data using the Soding 
parameterisation is shown in the figure, where contributions from the p° resonance, non-resonant 7r + 7r~ 
pair production and the interference term are displayed. Different ranges of the invariant mass M n +^- 
were examined in order to minimise the impact of non-resonant pair production and the interference 
term. The determined range is given by Eq. (|6]>, leading to an overall effect of less than 2 % in any 
kinematic bin of Q 2 , x Bj or pj,. Exclusive resonant p° and exclusive non-resonant pair contributions are 
not distinguished in the following. 

(iv) Coherent exclusive p° production on various nuclei of the target constitutes additional background. 
Its magnitude is estimated from the analysis of the shape of p\ distributions. In the kinematic region 
defined by Eqs (j6]>-([T0|, it amounts to ~ 12% for NH3 and ~ 8% for 6 LiD. No correction is applied 



for this residual background, which is justified by the observation that in this region the asymmetry 
is consistent within statistical uncertainty with that for events from the small-p^ region. The 
latter is defined by p\ < 0.05 (GeV/c) 2 for the NH 3 target and p\ < 0.10 (GeV/c) 2 for the 6 LiD target, 
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Fig. 7: The M n+7T ~ distribution with a fit using the Soding parameterisation is shown for NH3 data 
in the selected kinematic range 1 (GeV/c) 2 < Q 2 < 1.2(GeV/c) 2 , using the restrictive cut — 2.5GeV < 
-EWss < OGeV to suppress the semi-inclusive background. The thick solid line represents the result of the 
fit, while the dashed, dotted and dashed-dotted ones represent contributions from resonant p° production, 
non-resonant ir + ir~ production and the interference term respectively. 



where coherent production dominates. 



6 Extraction of the asymmetry 

The asymmetry is determined in bins of Q 2 , xbj or p 2 ^, while integrating over the two remaining vari- 
ables. For brevity the dependence on these variables is omitted in the following. The number of exclusive 
events, after subtraction of SIDIS background, can be expressed as a function of the angle 4> — 4>s in the 
following way 



N^-<p s )=Fnaa Q [l±f\P T \A^ VSI sin(0 - <fe)J , (12) 

where F is the muon flux, n the number of target nucleons, a{4> — 4>s) the product of acceptance and 
efficiency of the apparatus, gq the spin-averaged cross section, / the dilution factor, ±|Pt| the target 
polarisation, and the asymmetry A™^? is defined by Eq. ((5]). In this analysis, the asymmetry is 
extracted from a direct fit of the number of events in bins of 4> — <ps- For each of the two target cells 
(ncell = 1,2) and polarisation states (+, — ), 



^ncell = T / W 4,ncell(l ±Asin(0 ')) (13) 

J A -A/2 J 



gives the number of events in bin j of 4> — <^s, where j = 1 to m, and m = 12 the number of bins. Here 
4>' = 4> — 4>s, the central value of bin j is denoted by </> CJ - and A is the bin width. In the case of NH3 data, 
events from upstream and downstream cells of the three-cell target are combined. This leads to a system 
of 2 x 2 x m non-linear equations. The 'normalised acceptance' c^ ncell is the product of spin-averaged 
cross section, muon flux, number of target nucleons, acceptance and efficiency of the apparatus. The 
dependence on target dilution factor / and target polarisation Pt is absorbed into the 'raw asymmetry' 
A that is directly fitted to the data. 

A one-dimensional binned maximum likelihood method is used to solve the system of Am equations. 
Here the likelihood is constructed from Gaussian distributions in order to account for the non-Poissonian 
nature of the background subtracted data. In order to reduce the number of unknowns, the reasonable 
assumption is made that possible changes of acceptances in target cell before and after target polarisation 
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6 EXTRACTION OF THE ASYMMETRY 



reversal are the same for every bin j and can be described by a common constant C: 

c + c + 

C= J±J^2_ (14) 

Using this constraint, one can determine the 3m independent normalised acceptances c^ ncell , the constant 
C and the raw asymmetry A. 

In every bin XBj, Q 2 or p^, the asymmetry A^ff ^ s hs calculated as A & ™^ ^ s '= A/(f \Pt\) using 
the raw asymmetry A obtained from the fit. The dilution factor / is calculated on an event-by-event 
basis using the measured contributions of various atomic elements in the target and a parameterisation of 
the nuclear dependence of the spin-independent cross section for the studied reaction, pN — > pp°N, as 
explained in Ref. ll32ll . The ratios of this cross section per nucleon for a given nucleus to the cross section 
on the proton or deuteron are parameterised [33 ] over a wide Q 2 -range, using measurements on various 
nuclear targets. No dependence of nuclear effects on y or v is assumed, motivated by NMC results on 
exclusive p° production 11341 in a kinematic range similar to that of COMPASS. The Q 2 dependence of 
the dilution factor for incoherent exclusive p° production is shown in Fig. [8] As can be seen, the values of 
the dilution factor / for the NH 3 target vary from 0.27 at Q 2 = 1 (GeV/c) 2 to 0.18 at Q 2 = 10 (GeV/c) 2 , 
and correspondingly from 0.45 to 0.42 for the 6 LiD targefj] Radiative corrections are neglected in the 
present analysis, in particular in the calculation of /. They are expected to be small, mainly because of the 
exclusivity cuts (see Sec. [4]) that largely suppress the otherwise dominant external photon radiation 1351 . 
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Fig. 8: Dilution factor / for exclusive p° production as a function of Q 2 . The solid (dashed) line repre- 
sents the dilution factor for the NH3 ( 6 LiD) target. 



The correction to the proton asymmetry A V j, due to the polarisation of N nuclei in the ammonia 
target was estimated following the approach of Refs ||3"6l [37J. The correction is proportional to the 
measured asymmetry for the deuteron and approximately given by: 

. ,sin(0-</> s ) _ 1 / 1-v 1 a d ,sin(</>-0 s ) /in 

aa ut, p -3- ( ~3 j '6'^ UT ' d ■ ( 5) 

The factors account for the fraction of polarisable nitrogen nuclei in ammonia, the alignment of proton 
spin vs. 14 N spin, the ratio of 14 N to l H polarisations and the ratio of cross sections, 0"^ and a p , for 
exclusive p° production by scattering muons off unpolarised deuteron and proton targets, respectively. 
The estimated corrections are very small, typically about 0.1 %, and are neglected in the following. 

^ur estimates of / for the 6 LiD target are by 9 - 19 % higher than those of Ref. 1321 . The difference, which is due to 
an improvement of the treatment of the 6 Li nucleus, is significantly smaller than the total systematic uncertainty quoted in 
Ref. 03. 
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7 Systematic uncertainties 

In this section we describe tests performed to examine various sources of possible systematic uncertain- 
ties, namely: (a) bias of the applied estimator of the asymmetry, (b) data stability, (c) false asymmetries, 
(d) sensitivity to the method of background subtraction, (e) sensitivity to the Monte Carlo description 
used for the parameterisation of the background shape, (f) compatibility of results after background 
subtraction, and (g) uncertainties of target dilution factor and target polarisation value. As tests (a)-(c) 
are not sensitive to the background, they are performed using background non-corrected data. Possible 
systematic effects related to background subtraction are subject of tests (d)-(f). 

(a) In this analysis a one-dimensional binned maximum likelihood method is used (see Sec. [6]) rather than 
an extended unbinned maximum likelihood method IT381 l39ll because the latter one depends strongly on 
the quality of the Monte Carlo description for the SIDIS background in the kinematic region of our data. 
At present none of the existing generators satisfactorily describes the SIDIS background. The bin-by -bin 
comparison of the background non-corrected results for each target indicates good agreement between 
both estimators. However, for each of the three binnings used, in XBj, Q 2 and p\, the mean asymmetry 
values for the deuteron from the binned method is observed to be slightly smaller than that from the 
unbinned method, by several percent of the statistical uncertainty. Additionally, systematic differences 
of similar size are seen in the deuteron data between the mean values of the asymmetries evaluated in the 
three binnings. Both these effects result in a systematic uncertainty of 0. 10 a stat . As systematic uncer- 
tainties (a)-(f) are evaluated using data, here and in the following a systematic uncertainty is expressed 
in terms of a stat , which is the statistical error of the background corrected asymmetry measured in a 
given kinematic bin. For the proton data one observes good agreement between different estimators also 
at the level of mean asymmetries. The distribution of differences between mean asymmetries obtained 
from two estimators is centred at zero and its RMS value is used to estimate a systematic uncertainty of 
«0. 12 o- stai . 

(b) Subsets of data, each containing two periods consecutive in time with opposite target polarisations, 
are compared to test the stability of data taking. For each of the 18 (3) subsets formed from proton 
(deuteron) data, the asymmetry is determined in every kinematic bin of x Bj, Q 2 and p\. All asymmetries 
are found to be compatible within statistical uncertainties. 

(c) In order to investigate possible false asymmetries, the target is artificially divided into four cells 
of 30 cm length each, distributed contiguously along the beam direction. This allows the evaluation 
of two independent false asymmetries using cells with the same spin orientation. Determining these 
false asymmetries in the exclusive region leads to statistical fluctuations similar to those of the physics 
asymmetry. In order to increase the statistical significance, the false asymmetries are studied in an 
extended range, — 10GeV< -E m i ss < 20GeV. The resulting false asymmetries are found to be consistent 
with zero within statistical uncertainties. They are used to estimate upper bounds for the corresponding 
systematic uncertainty, namely 0.15<7 s * a * (0.49cr stot ) for the proton (deuteron) data. 

(d) In order to estimate the sensitivity of the extracted azimuthal asymmetry to the method of back- 
ground subtraction, an alternative method is applied using measured 4> — 4>s distributions of the SIDIS 
background. Their shapes are determined in the range 7 GeV < E m { ss < 20 GeV, where exclusive and 
diffractive-dissociation events can be neglected against the SIDIS background. They are measured for 
each kinematic bin in XBj, Q 2 , and p\ respectively, and for each target cell and polarisation state. They 
are rescaled in such a way that in the exclusive region (see Eq. ([7])) the total number of events for the 
rescaled distribution is equal to the number of background events obtained from the signal plus back- 
ground fit to the -E m i ss spectra. The rescaled spectra are subtracted from the corresponding <j) — 4>s 
distributions for the data and the asymmetries are extracted from the resulting background-corrected 
distributions as described in Sec. [6] In this approach one relies on the assumption that the background 
asymmetry does not depend on i?miss- It is supported by studies of the asymmetry as a function of -E m i ss hi 
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7 SYSTEMATIC UNCERTAINTIES 



Table 3: Estimates of systematic uncertainties of A^f ^for proton and deuteron data. No value is 
quoted when the systematic uncertainty is negligibly small, below 0.01. 

source of systematic uncertainty a sys /o stat a sys /a stat 

proton data deuteron data 

(a) ID binned estimator 0.12 0.10 

(b) data stability — — 

(c) false asymmetries 0.15 0.49 

(d) method of background subtraction — — 

(e) MC dependence — 0.04 

(f) compatibility after background subtraction in xsj, Q and p\ 0.15 — 

total 0.25 0.50 



the range — 2.5GeV < EWss < 20GeV where no such dependence was found within statistical precision. 
Both the default and alternative methods are used to extract A*™^ A point-by-point comparison in 
the three kinematic binnings indicates for all data sets very good agreement between the two methods, 
within statistical uncertainties. 

(e) In order to study the sensitivity of the SIDIS background subtraction to the Monte Carlo description, 
the effect of using unweighted or weighted LEPTO samples for the background parameterisation is 
investigated. Despite the different shapes of -E m i ss distributions in the two cases, differences between 
background-corrected asymmetries are very small as expected, because a bad parameterisation of the 
background does not introduce a background asymmetry. Additionally, for the 2007 set-up a second large 
Monte Carlo background sample was generated using PYTHIA with default values of parameters [40]. 
The background shapes obtained with weighted LEPTO and weighted PYTHIA are similar but in general 
the latter results in about 10% less background. In most of the kinematic bins the asymmetries are 
the same. The systematic uncertainty due to the Monte Carlo description of the SIDIS background is 
negligible for the proton data and is estimated to be about 0.04<r s ' a * for the deuteron data. 

(f) An important consistency test is the comparison of the mean asymmetry values evaluated in bins of 
XBj, Q 2 and pj,, to check if the assumption of Eq. (hfy holds after background subtraction. For the 
proton data, the mean values of the asymmetries evaluated in bins of xsj and Q 2 are compatible, and the 
one evaluated in bins of p\ agrees within about 030a stat ; half of this difference is taken into account as 
systematic uncertainty due to background subtraction. For the deuteron data, the three mean asymmetry 
values are in reasonable agreement. The small differences observed are similar as in the case of the 
background non-corrected asymmetries and hence not introduced by background subtraction. 

A summary of the systematic uncertainties (a)-(f) is given in Table [3] For each kinematic bin i, the total 
systematic uncertainty of the measured asymmetry A*™j? ^ s 'i s obtained as quadratic superposition of 
the sources (a) to (f). It equals a* ya = 0.25 af tat for the proton and af 3 = 0.5 af tat for the deuteron data, 

where af tat is the statistical uncertainty of the measured asymmetry Afjtf in bin i. 

Not listed in the table are the scale uncertainties due to the relative uncertainties of dilution factor and 
target polarisation. For the NH3 target these are 2% and 3%, which leads by quadratic superposition to 
a scale uncertainty of 0.036. The analogous values for the 6 LiD target are 2% and 5% , respectively, 
resulting in a scale uncertainty of 0.054. Note that these scale uncertainties are common for all measured 
asymmetries for a given target. 
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8 Results and comparison to model predictions 

The transverse target spin asymmetries A s ^jf ^ measured on proton and deuteron are shown in Fig. [9] 
as a function of xbj, Q 2 or pj,, upon integrating over the two other variables. For both targets the 
asymmetries are found to be small and consistent with zero within statistical uncertainties. Note that 
this is the first measurement of A s ^f ^for transversely polarised deuterons. The numerical values for 
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Fig. 9: Transverse target spin asymmetries A s ™^ ™ measured on proton (upper) and deuteron (lower) 
as a function of xsj, Q 2 and p\. Error bars show statistical uncertainties, while the systematic ones are 
represented by grey bands at the bottom. The curves show the predictions of the GPD model [41 J using 
the set of parameters called 'variant 1'. They are calculated at W = 8.1 GeV/c 2 and p\ = 0.2(GeV/c) 
for the left and middle panels, and at W = 8.1 GeV/c 2 and Q 2 = 2.2 (GeV/c) 2 for the right panels. The 
theoretical error bands reflect uncertainties of GPD parameterisations. 



j^m[<p vs)^ presented in Table |4| for each XBj, Q 2 and p\ bin, together with statistical and system- 
atic uncertainties. Also, average values of kinematic variables for each bin are given. Averaged over 
the COMPASS kinematic region, the values of A^~^\xe -0.002 ±0.010(stat) ±0.003(sys) for the 
proton and 0.02 ±0.03(stat) ±0.02(sys) for the deuteron. 

The results of a similar measurement of the asymmetry A™^? for p° production on the proton target 
by the HERMES experiment lfl9l are also consistent with zero within total experimental uncertainties. 
The separate asymmetries for longitudinally and transversely polarised p° mesons were found by HER- 
MES [ 20] to be consistent with zero as well. 

Theoretical predictions for A*™j? for p° are given by the GPD model of Goloskokov and Kroll Il4~lll . 
In this model, electroproduction of a light vector meson V at small x Bj is analysed in the handbag ap- 
proach, in which the amplitude of the process is a convolution of GPDs with amplitudes for the partonic 
subprocesses 7*g^ — > Vq? and j*g — > Vg. The partonic subprocess amplitudes, which comprise corre- 
sponding hard scattering kernels and meson DAs, are calculated in the modified perturbative approach 
where the transverse momenta of quarks and antiquarks forming the vector meson are retained and Su- 
dakov suppressions are taken into account. The model gives predictions for contributions from both 



18 



8 RESULTS AND COMPARISON TO MODEL PREDICTIONS 



Table 4: The transverse target spin asymmetries A v ^ measured on proton and deuteron in bins 
of Q 2 , XBj and pj,. The systematic uncertainties are obtained using the values given in Table [3] In 
addition, a scale uncertainty of 3.6 % (5.4 %) accounts for uncertainties in the determination of the target 
polarisation and target dilution factor for proton (deuteron) data. 



proton data 





(Q 2 ) 


(GeV/c) 2 


( x Bj) 


(P 2 T ) 


(GeV/c) 2 


A sin((j>-(f>s) 1 _ s tat 1 „sys 
A UT,p ±a ±G 


Q 2 bin (GeV/c) 2 


1.0-1.2 




1.1 


0.019 




0.18 


-0.019 ±0.019 ±0.005 


1.2-1.6 




1.4 


0.025 




0.18 


-0.003 ±0.018 ±0.004 


1.6-2.4 




1.9 


0.035 




0.18 


-0.001 ±0.020 ±0.005 


2.4-10.0 




4.0 


0.076 




0.19 


0.027 ±0.026 ±0.007 


xbj bin 


0.003-0.02 




1.4 


0.013 




0.17 


-0.010 ±0.015 ±0.004 


0.02-0.03 




1.6 


0.025 




0.18 


0.018 ±0.020 ±0.005 


0.03-0.05 




1.9 


0.038 




0.18 


-0.015 ±0.020 ±0.005 


0.05-0.30 




3.8 


0.088 




0.19 


0.015 ±0.031 ±0.008 


p 2 T bin (GeV/c) 2 


0.05-0.10 




2.1 


0.037 




0.07 


0.000 ±0.016 ±0.004 


0.10-0.15 




2.1 


0.039 




0.12 


0.017 ±0.023 ±0.006 


0.15-0.25 




2.2 


0.040 




0.20 


0.001 ±0.022 ±0.005 


0.25-0.35 




2.2 


0.042 




0.30 


-0.039 ±0.031 ±0.008 


0.35-0.50 




2.3 


0.043 




0.42 


0.012 ±0.036 ±0.009 








deuteron data 








(Q 2 ) 


(GeV/c) 2 


( x Bj) 


(p 2 t) 


(GeV/c) 2 


4 sin(0-</> s ) , stat , sys 
■^UT,d ^° ^° 


Q 2 bin (GeV/c) 2 


1.0- 1.2 




1.1 


0.018 




0.23 


0.09 ±0.05 ±0.03 


1.2-1.6 




1.4 


0.023 




0.23 


0.00 ±0.05 ±0.02 


1.6-2.4 




1.9 


0.031 




0.23 


-0.03 ±0.06 ±0.03 


2.4-10.0 




3.9 


0.059 




0.24 


0.00 ±0.07 ±0.03 


xbj bin 


0.003-0.02 




1.4 


0.013 




0.23 


0.03 ±0.04 ±0.02 


0.02-0.03 




1.6 


0.025 




0.23 


0.01 ±0.06 ±0.03 


0.03-0.05 




2.0 


0.038 




0.23 


-0.02 ±0.06 ±0.03 


0.05-0.30 




3.9 


0.078 




0.24 


0.04 ±0.09 ±0.05 


p 2 T bin (GeV/c) 2 


0.10-0.15 




1.9 


0.031 




0.12 


0.02 ±0.05 ±0.02 


0.15-0.25 




2.0 


0.031 




0.19 


0.04 ±0.05 ±0.02 


0.25-0.35 




2.0 


0.032 




0.30 


-0.01 ±0.07 ±0.03 


0.35-0.50 




2.1 


0.033 




0.42 


-0.03 ±0.08 ±0.04' 
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longitudinal and transverse virtual photons. 

The predicted value of the proton asymmetry averaged over the COMPASS kinematic region is about 
—0.02 and correspondingly about —0.03 for HERMES. The comparison of the COMPASS results as 
a function of xbj, Q 2 and p\* to the predictions of the model, shown in Fig. |9j indicates reasonable 
agreement. The curves were obtained using the default version of the model ('variant 1'), with only 
contributions from valence quark GPDs E u and E d . The indicated theoretical error bands reflect uncer- 
tainties in the GPD parameterisations. 

In order to investigate the role of gluons and sea quarks the authors of the model consider two extreme 
cases for non-zero GPDs E 9 and E sea . They use either positive or negative E sea that saturates positivity 
bounds, and E 9 that is constrained by a sum-rule for the second moments of GPDs E of quarks and 
gluons fiTll . Including GPDs E 9 and E sea has a very small effect on the predicted values of A s ™j? 
resulting in differences with respect to the default version which are significantly smaller than the theo- 
retical uncertainties shown in Fig. [9] 

The sensitivity of A s ™j? to the light quark GPDs, E u and E d , is different for the two targets. For 
the proton these GPDs enter the amplitude as the sum 2/3 E u + 1/3 E d . For incoherent production on 
the nucleons of the deuteron, assuming isospin invariance and neglecting nuclear effects, they effectively 
contribute as E u + E d . In both cases, a small value of the asymmetry for p° is expected as E u and E d 
are similar in magnitude but of opposite sign. 

In conclusion, the transverse target spin asymmetry A s ^jf~^for hard exclusive p° meson production 
was measured at COMPASS on the proton and, for the first time, on the deuteron. The values of the 
asymmetry for both targets are small and compatible with zero in a broad kinematic range. They are 
compatible with the predictions of the GPD model of Ref. ATI . The COMPASS proton results are in 
good agreement with those obtained at HERMES, while they are more precise by a factor of about 3 and 
cover a larger kinematic domain. 

We acknowledge the support of the CERN management and staff, as well as the skills and efforts of the 
technicians of the collaborating institutes. We also thank S.V. Goloskokov and P. Kroll for discussions 
of the results. 
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